High-throughput mRNA sequencing of human placenta shows sex differences across gestation.

INTRODUCTION: Fetal sex affects fetal and maternal health outcomes in pregnancy, but this connection remains poorly understood. As the placenta is the route of fetomaternal communication and derives from the fetal genome, placental gene expression sex differences may explain these outcomes. OBJECTIVES: We utilized next generation sequencing to study the normal human placenta in both sexes in first and third trimester to generate a normative transcriptome based on sex and gestation. STUDY DESIGN: We analyzed 124 first trimester (T1, 59 female and 65 male) and 43 third trimester (T3, 18 female and 25 male) samples for sex differences within each trimester and sex-specific gestational differences. RESULTS: Placenta shows more significant sexual dimorphism in T1, with 94 T1 and 26 T3 differentially expressed genes (DEGs). The sex chromosomes contributed 60.6% of DEGs in T1 and 80.8% of DEGs in T3, excluding X/Y pseudoautosomal regions. There were 6 DEGs from the pseudoautosomal regions, only significant in T1 and all upregulated in males. The distribution of DEGs on the X chromosome suggests genes on Xp (the short arm) may be particularly important in placental sex differences. Dosage compensation analysis of X/Y homolog genes shows expression is primarily contributed by the X chromosome. In sex-specific analyses of first versus third trimester, there were 2815 DEGs common to both sexes upregulated in T1, and 3263 common DEGs upregulated in T3. There were 7 female-exclusive DEGs upregulated in T1, 15 female-exclusive DEGs upregulated in T3, 10 male-exclusive DEGs upregulated in T1, and 20 male-exclusive DEGs upregulated in T3. DISCUSSION: This is the largest cohort of placentas across gestation from healthy pregnancies defining the normative sex dimorphic gene expression and sex common, sex specific and sex exclusive gene expression across gestation. The first trimester has the most sexually dimorphic transcripts, and the majority were upregulated in females compared to males in both trimesters. The short arm of the X chromosome and the pseudoautosomal region is particularly critical in defining sex differences in the first trimester placenta. As pregnancy is a dynamic state, sex specific DEGs across gestation may contribute to sex dimorphic changes in overall outcomes.

Fallopian tube single cell analysis reveals myeloid cell alterations in high-grade serous ovarian cancer.

Most high-grade serous ovarian cancers (HGSCs) likely initiate from fallopian tube (FT) epithelia. While epithelial subtypes have been characterized using single-cell RNA-sequencing (scRNA-Seq), heterogeneity of other compartments and their involvement in tumor progression are poorly defined. Integrated analysis of human FT scRNA-Seq and HGSC-related tissues, including tumors, revealed greater immune and stromal transcriptional diversity than previously reported. We identified abundant monocytes in FTs across two independent cohorts. The ratio of macrophages to monocytes is similar between benign FTs, ovaries, and adjacent normal tissues but significantly greater in tumors. FT-defined monocyte and macrophage signatures, cell-cell communication, and gene set enrichment analyses identified monocyte- and macrophage-specific interactions and functional pathways in paired tumors and adjacent normal tissues. Further reanalysis of HGSC scRNA-Seq identified monocyte and macrophage subsets associated with neoadjuvant chemotherapy. Taken together, our work provides data that an altered FT myeloid cell composition could inform the discovery of early detection markers for HGSC.

MYBL2 drives prostate cancer plasticity and identifies CDK2 as a therapeutic vulnerability in RB1-loss and neuroendocrine prostate cancer.

Phenotypic plasticity is a recognized mechanism driving therapeutic resistance in prostate cancer (PCa) patients. While underlying molecular causations driving phenotypic plasticity have been identified, therapeutic success is yet to be achieved. To identify putative master regulator transcription factors (MR-TF) driving phenotypic plasticity in PCa, this work utilized a multiomic approach using genetically engineered mouse models of prostate cancer combined with patient data to identify MYBL2 as a significantly enriched transcription factor in PCa exhibiting phenotypic plasticity. Genetic inhibition of Mybl2 using independent murine PCa cell lines representing phenotypic plasticity demonstrated Mybl2 loss significantly decreased in vivo growth as well as cell fitness and repressed gene expression signatures involved in pluripotency and stemness. Because MYBL2 is currently not druggable, a MYBL2 gene signature was employed to identify cyclin-dependent kinase-2 (CDK2) as a potential therapeutic target. CDK2 inhibition phenocopied genetic loss of Mybl2 and significantly decreased in vivo tumor growth associated with enrichment of DNA damage. Together, this work demonstrates MYBL2 as an important MR-TF driving phenotypic plasticity in PCa. Further, high MYBL2 activity identifies PCa that would be responsive to CDK2 inhibition. Significance: PCa that escapes therapy targeting the androgen receptor signaling pathways via phenotypic plasticity are currently untreatable. Our study identifies MYBL2 as a MR-TF in phenotypic plastic PCa and implicates CDK2 inhibition as novel therapeutic target for this most lethal subtype of PCa.

High-throughput mRNA-seq atlas of human placenta shows vast transcriptome remodeling from first to third trimester†.

The placenta, composed of chorionic villi, changes dramatically across gestation. Understanding differences in ongoing pregnancies are essential to identify the role of chorionic villi at specific times in gestation and develop biomarkers and prognostic indicators of maternal-fetal health. The normative mRNA profile is established using next-generation sequencing of 124 first trimester and 43 third trimester human placentas from ongoing healthy pregnancies. Stably expressed genes (SEGs) not different between trimesters and with low variability are identified. Differential expression analysis of first versus third trimester adjusted for fetal sex is performed, followed by a subanalysis with 23 matched pregnancies to control for subject variability using the same genetic and environmental background. Placenta expresses 14,979 polyadenylated genes above sequencing noise (transcripts per million > 0.66), with 10.7% SEGs across gestation. Differentially expressed genes (DEGs) account for 86.7% of genes in the full cohort [false discovery rate (FDR) < 0.05]. Fold changes highly correlate between the full cohort and subanalysis (Pearson = 0.98). At stricter thresholds (FDR < 0.001, fold change > 1.5), there remains 50.1% DEGs (3353 upregulated in first and 4155 upregulated in third trimester). This is the largest mRNA atlas of healthy human placenta across gestation, controlling for genetic and environmental factors, demonstrating substantial changes from first to third trimester in chorionic villi. Specific differences and SEGs may be used to understand the specific role of the chorionic villi throughout gestation and develop first trimester biomarkers of placental health that transpire across gestation, which can be used for future development of biomarkers for maternal-fetal health.

α-Ketoglutarate-Dependent KDM6 Histone Demethylases and Interferon-Stimulated Gene Expression in Lupus.

OBJECTIVE: We aimed to investigate the hypothesis that interferon (IFN)-stimulated gene (ISG) expression in systemic lupus erythematosus (SLE) monocytes is linked to changes in metabolic reprogramming and epigenetic regulation of ISG expression. METHODS: Monocytes from healthy volunteers and patients with SLE at baseline or following IFNα treatment were analyzed by extracellular flux analysis, proteomics, metabolomics, chromatin immunoprecipitation, and gene expression. The histone demethylases KDM6A/B were inhibited using glycogen synthase kinase J4 (GSK-J4). GSK-J4 was tested in pristane and resiquimod (R848) models of IFN-driven SLE. RESULTS: SLE monocytes had enhanced rates of glycolysis and oxidative phosphorylation compared to healthy control monocytes, as well as increased levels of isocitrate dehydrogenase and its product, α-ketoglutarate (α-KG). Because α-KG is a required cofactor for histone demethylases KDM6A and KDM6B, we hypothesized that IFNα may be driving "trained immune" responses through altering histone methylation. IFNα priming (day 1) resulted in a sustained increase in the expression of ISGs in primed cells (day 5) and enhanced expression on restimulation with IFNα. Importantly, decreased H3K27 trimethylation was observed at the promoters of ISGs following IFNα priming. Finally, GSK-J4 (KDM6A/B inhibitor) resulted in decreased ISG expression in SLE patient monocytes, as well as reduced autoantibody production, ISG expression, and kidney pathology in R848-treated BALB/c mice. CONCLUSION: Our study suggests long-term IFNα exposure alters the epigenetic regulation of ISG expression in SLE monocytes via changes in immunometabolism, a mechanism reflecting trained immunity to type I IFN. Importantly, it opens the possibility that targeting histone-modifying enzymes, such as KDM6A/B, may reduce IFN responses in SLE.

The tumor multi-omic landscape of endometrial cancers developed on a germline genetic background of adiposity.

High body mass index (BMI) is a causal risk factor for endometrial cancer but the tumor molecular mechanisms affected by adiposity and their therapeutic relevance remain poorly understood. Here we characterize the tumor multi-omic landscape of endometrial cancers that have developed on a background of lifelong germline genetic exposure to elevated BMI. We built a polygenic score (PGS) for BMI in women using data on independent, genome-wide significant variants associated with adult BMI in 434,794 women. We performed germline (blood) genotype quality control and imputation on data from 354 endometrial cancer cases from The Cancer Genome Atlas (TCGA). We assigned each case in this TCGA cohort their genetically predicted life-course BMI based on the BMI PGS. Multivariable generalized linear models adjusted for age, stage, microsatellite status and genetic principal components were used to test for associations between the BMI germline PGS and endometrial cancer tumor genome-wide genomic, transcriptomic, proteomic, epigenomic and immune traits in TCGA. High BMI germline PGS was associated with (i) upregulated tumor gene expression in the IL6-JAK-STAT3 pathway (FDR=4.2×10-7); (ii) increased estimated intra-tumor activated mast cell infiltration (FDR=0.008); (iii) increased single base substitution (SBS) mutational signatures 1 (FDR=0.03) and 5 (FDR=0.09) and decreased SBS13 (FDR=0.09), implicating age-related and APOBEC mutagenesis, respectively; and (iv) decreased tumor EGFR protein expression (FDR=0.07). Alterations in IL6-JAK-STAT3 signaling gene and EGFR protein expression were, in turn, significantly associated with both overall survival and progression-free interval. Thus, we integrated germline and somatic data using a novel study design to identify associations between genetically predicted lifelong exposure to higher BMI and potentially actionable endometrial cancer tumor molecular features. These associations inform our understanding of how high BMI may influence the development and progression of this cancer, impacting endometrial tumor biology and clinical outcomes.

A biallelic multiple nucleotide length polymorphism explains functional causality at 5p15.33 prostate cancer risk locus.

To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

High-throughput mRNA-seq atlas of human placenta shows vast transcriptome remodeling from first to third trimester.

Background: The placenta, composed of chorionic villi, changes dramatically across gestation. Understanding differences in ongoing pregnancies are essential to identify the role of chorionic villi at specific times in gestation and develop biomarkers and prognostic indicators of maternal- fetal health. Methods: The normative mRNA profile is established using next-generation sequencing of 124 first trimester and 43 third trimester human placentas from ongoing healthy pregnancies. Stably expressed genes not different between trimesters and with low variability are identified. Differential expression analysis of first versus third trimester adjusted for fetal sex is performed, followed by a subanalysis with 23 matched pregnancies to control for subject variability using the same genetic and environmental background. Results: Placenta expresses 14,979 mRNAs above sequencing noise (TPM>0.66), with 1,545 stably expressed genes across gestation. Differentially expressed genes account for 86.7% of genes in the full cohort (FDR<0.05). Fold changes highly correlate between the full cohort and subanalysis (Pearson = 0.98). At stricter thresholds (FDR<0.001, fold change>1.5), there are 6,941 differentially expressed protein coding genes (3,206 upregulated in first and 3,735 upregulated in third trimester). Conclusion: This is the largest mRNA atlas of healthy human placenta across gestation, controlling for genetic and environmental factors, demonstrating substantial changes from first to third trimester in chorionic villi. Specific differences and stably expressed genes may be used to understand the specific role of the chorionic villi throughout gestation and develop first trimester biomarkers of placental health that transpire across gestation, which can be used for future development of biomarkers in maternal-fetal disease.

Rewiring of master transcription factor cistromes during high-grade serous ovarian cancer development.

The transcription factors MECOM, PAX8, SOX17 and WT1 are candidate master regulators of high-grade serous 'ovarian' cancer (HGSC), yet their cooperative role in the hypothesized tissue of origin, the fallopian tube secretory epithelium (FTSEC) is unknown. We generated 26 epigenome (CUT&TAG, CUT&RUN, ATAC-seq and HiC) data sets and 24 profiles of RNA-seq transcription factor knock-down followed by RNA sequencing in FTSEC and HGSC models to define binding sites and gene sets regulated by these factors in cis and trans. This revealed that MECOM, PAX8, SOX17 and WT1 are lineage-enriched, super-enhancer associated master regulators whose cooperative DNA-binding patterns and target genes are re-wired during tumor development. All four TFs were indispensable for HGSC clonogenicity and survival but only depletion of PAX8 and WT1 impaired FTSEC cell survival. These four TFs were pharmacologically inhibited by transcriptional inhibitors only in HGSCs but not in FTSECs. Collectively, our data highlights that tumor-specific epigenetic remodeling is tightly related to MECOM, PAX8, SOX17 and WT1 activity and these transcription factors are targetable in a tumor-specific manner through transcriptional inhibitors.

KRAS mutations and endometriosis burden of disease.

The clinical phenotype of somatic mutations in endometriosis is unknown. The objective was to determine whether somatic KRAS mutations were associated with greater disease burden in endometriosis (i.e. more severe subtypes and higher stage). This prospective longitudinal cohort study included 122 subjects undergoing endometriosis surgery at a tertiary referral center between 2013 and 2017, with 5-9 years of follow-up. Somatic activating KRAS codon 12 mutations were detected in endometriosis lesions using droplet digital PCR. KRAS mutation status for each subject was coded as present (KRAS mutation in at least one endometriosis sample in a subject) or absent. Standardized clinical phenotyping for each subject was carried out via linkage to a prospective registry. Primary outcome was anatomic disease burden, based on distribution of subtypes (deep infiltrating endometriosis, ovarian endometrioma, and superficial peritoneal endometriosis) and surgical staging (Stages I-IV). Secondary outcomes were markers of surgical difficulty, demographics, pain scores, and risk of re-operation. KRAS mutation presence was higher in subjects with deep infiltrating endometriosis or endometrioma lesions only (57.9%; 11/19) and subjects with mixed subtypes (60.6%; 40/66), compared with those with superficial endometriosis only (35.1%; 13/37) (p = 0.04). KRAS mutation was present in 27.6% (8/29) of Stage I cases, in comparison to 65.0% (13/20) of Stage II, 63.0% (17/27) of Stage III, and 58.1% (25/43) of Stage IV cases (p = 0.02). KRAS mutation was also associated with greater surgical difficulty (ureterolysis) (relative risk [RR] = 1.47, 95% CI: 1.02-2.11) and non-Caucasian ethnicity (RR = 0.64, 95% CI: 0.47-0.89). Pain severities did not differ based on KRAS mutation status, at either baseline or follow-up. Re-operation rates were low overall, occurring in 17.2% with KRAS mutation compared with 10.3% without (RR = 1.66, 95% CI: 0.66-4.21). In conclusion, KRAS mutations were associated with greater anatomic severity of endometriosis, resulting in increased surgical difficulty. Somatic cancer-driver mutations may inform a future molecular classification of endometriosis.

Epigenomic charting and functional annotation of risk loci in renal cell carcinoma.

While the mutational and transcriptional landscapes of renal cell carcinoma (RCC) are well-known, the epigenome is poorly understood. We characterize the epigenome of clear cell (ccRCC), papillary (pRCC), and chromophobe RCC (chRCC) by using ChIP-seq, ATAC-Seq, RNA-seq, and SNP arrays. We integrate 153 individual data sets from 42 patients and nominate 50 histology-specific master transcription factors (MTF) to define RCC histologic subtypes, including EPAS1 and ETS-1 in ccRCC, HNF1B in pRCC, and FOXI1 in chRCC. We confirm histology-specific MTFs via immunohistochemistry including a ccRCC-specific TF, BHLHE41. FOXI1 overexpression with knock-down of EPAS1 in the 786-O ccRCC cell line induces transcriptional upregulation of chRCC-specific genes, TFCP2L1, ATP6V0D2, KIT, and INSRR, implicating FOXI1 as a MTF for chRCC. Integrating RCC GWAS risk SNPs with H3K27ac ChIP-seq and ATAC-seq data reveals that risk-variants are significantly enriched in allelically-imbalanced peaks. This epigenomic atlas in primary human samples provides a resource for future investigation.

Single-cell transcriptomic analysis of endometriosis.

Endometriosis is a common condition in women that causes chronic pain and infertility and is associated with an elevated risk of ovarian cancer. We profiled transcriptomes of >370,000 individual cells from endometriomas (n = 8), endometriosis (n = 28), eutopic endometrium (n = 10), unaffected ovary (n = 4) and endometriosis-free peritoneum (n = 4), generating a cellular atlas of endometrial-type epithelial cells, stromal cells and microenvironmental cell populations across tissue sites. Cellular and molecular signatures of endometrial-type epithelium and stroma differed across tissue types, suggesting a role for cellular restructuring and transcriptional reprogramming in the disease. Epithelium, stroma and proximal mesothelial cells of endometriomas showed dysregulation of pro-inflammatory pathways and upregulation of complement proteins. Somatic ARID1A mutation in epithelial cells was associated with upregulation of pro-angiogenic and pro-lymphangiogenic factors and remodeling of the endothelial cell compartment, with enrichment of lymphatic endothelial cells. Finally, signatures of ciliated epithelial cells were enriched in ovarian cancers, reinforcing epidemiologic associations between these two diseases.

circCsnk1g3- and circAnkib1-regulated interferon responses in sarcoma promote tumorigenesis by shaping the immune microenvironment.

Exonic circular RNAs (circRNAs) produce predominantly non-coding RNA species that have been recently profiled in many tumors. However, their functional contribution to cancer progression is still poorly understood. Here, we identify the circRNAs expressed in soft tissue sarcoma cells and explore how the circRNAs regulate sarcoma growth in vivo. We show that circCsnk1g3 and circAnkib1 promote tumor growth by shaping a pro-tumorigenic microenvironment, possibly due to their capabilities to regulate tumor-promoting elements extrinsic to the tumor cells. Accordingly, circCsnk1g3 and circAnkib1 can control the expression of interferon-related genes and pro-inflammatory factors in the sarcoma cells, thus directing immune cell recruitment into the tumor mass, and hence their activation. Mechanistically, circRNAs may repress pro-inflammatory elements by buffering activation of the pathways mediated by RIG-I, the cytosolic viral RNA sensor. The current findings suggest that the targeting of specific circRNAs could augment the efficacy of tumor and immune response to mainstay therapies.

Molecular Subclasses of Clear Cell Ovarian Carcinoma and Their Impact on Disease Behavior and Outcomes.

PURPOSE: To identify molecular subclasses of clear cell ovarian carcinoma (CCOC) and assess their impact on clinical presentation and outcomes. EXPERIMENTAL DESIGN: We profiled 421 primary CCOCs that passed quality control using a targeted deep sequencing panel of 163 putative CCOC driver genes and whole transcriptome sequencing of 211 of these tumors. Molecularly defined subgroups were identified and tested for association with clinical characteristics and overall survival. RESULTS: We detected a putative somatic driver mutation in at least one candidate gene in 95% (401/421) of CCOC tumors including ARID1A (in 49% of tumors), PIK3CA (49%), TERT (20%), and TP53 (16%). Clustering of cancer driver mutations and RNA expression converged upon two distinct subclasses of CCOC. The first was dominated by ARID1A-mutated tumors with enriched expression of canonical CCOC genes and markers of platinum resistance; the second was largely comprised of tumors with TP53 mutations and enriched for the expression of genes involved in extracellular matrix organization and mesenchymal differentiation. Compared with the ARID1A-mutated group, women with TP53-mutated tumors were more likely to have advanced-stage disease, no antecedent history of endometriosis, and poorer survival, driven by their advanced stage at presentation. In women with ARID1A-mutated tumors, there was a trend toward a lower rate of response to first-line platinum-based therapy. CONCLUSIONS: Our study suggests that CCOC consists of two distinct molecular subclasses with distinct clinical presentation and outcomes, with potential relevance to both traditional and experimental therapy responsiveness. See related commentary by Lheureux, p. 4838.

chromMAGMA: regulatory element-centric interrogation of risk variants.

Candidate causal risk variants from genome-wide association studies reside almost exclusively in noncoding regions of the genome and innovative approaches are necessary to understand their biological function. Multi-marker analysis of genomic annotation (MAGMA) is a widely used program that nominates candidate risk genes by mapping single-nucleotide polymorphism summary statistics from genome-wide association studies to gene bodies. We augmented MAGMA to create chromatin-MAGMA (chromMAGMA), a method to nominate candidate risk genes based on the presence of risk variants within noncoding regulatory elements (REs). We applied chromMAGMA to a genetic susceptibility dataset for epithelial ovarian cancer (EOC), a rare gynecologic malignancy characterized by high mortality. This identified 155 unique candidate EOC risk genes across five EOC histotypes; 83% (105/127) of high-grade serous ovarian cancer risk genes had not previously been implicated in this EOC histotype. Risk genes nominated by chromMAGMA converged on mRNA splicing and transcriptional dysregulation pathways. chromMAGMA is a pipeline that nominates candidate risk genes through a gene regulation-focused approach and helps interpret the biological mechanism of noncoding risk variants for complex diseases.

DNA methylation and transcriptomic features are preserved throughout disease recurrence and chemoresistance in high grade serous ovarian cancers.

BACKGROUND: Little is known about the role of global DNA methylation in recurrence and chemoresistance of high grade serous ovarian cancer (HGSOC). METHODS: We performed whole genome bisulfite sequencing and transcriptome sequencing in 62 primary and recurrent tumors from 28 patients with stage III/IV HGSOC, of which 11 patients carried germline, pathogenic BRCA1 and/or BRCA2 mutations. RESULTS: Landscapes of genome-wide methylation (on average 24.2 million CpGs per tumor) and transcriptomes in primary and recurrent tumors showed extensive heterogeneity between patients but were highly preserved in tumors from the same patient. We identified significant differences in the burden of differentially methylated regions (DMRs) in tumors from BRCA1/2 compared to non-BRCA1/2 carriers (mean 659 DMRs and 388 DMRs in paired comparisons respectively). We identified overexpression of immune pathways in BRCA1/2 carriers compared to non-carriers, implicating an increased immune response in improved survival (P = 0.006) in these BRCA1/2 carriers. CONCLUSION: These findings indicate methylome and gene expression programs established in the primary tumor are conserved throughout disease progression, even after extensive chemotherapy treatment, and that changes in methylation and gene expression are unlikely to serve as drivers for chemoresistance in HGSOC.

A multi-level investigation of the genetic relationship between endometriosis and ovarian cancer histotypes.

Endometriosis is associated with increased risk of epithelial ovarian cancers (EOCs). Using data from large endometriosis and EOC genome-wide association meta-analyses, we estimate the genetic correlation and evaluate the causal relationship between genetic liability to endometriosis and EOC histotypes, and identify shared susceptibility loci. We estimate a significant genetic correlation (rg) between endometriosis and clear cell (rg = 0.71), endometrioid (rg = 0.48), and high-grade serous (rg = 0.19) ovarian cancer, associations supported by Mendelian randomization analyses. Bivariate meta-analysis identified 28 loci associated with both endometriosis and EOC, including 19 with evidence for a shared underlying association signal. Differences in the shared risk suggest different underlying pathways may contribute to the relationship between endometriosis and the different histotypes. Functional annotation using transcriptomic and epigenomic profiles of relevant tissues/cells highlights several target genes. This comprehensive analysis reveals profound genetic overlap between endometriosis and EOC histotypes with valuable genomic targets for understanding the biological mechanisms linking the diseases.

The transcription factor PAX8 promotes angiogenesis in ovarian cancer through interaction with SOX17.

PAX8 is a master transcription factor that is essential during embryogenesis and promotes neoplastic growth. It is expressed by the secretory cells lining the female reproductive tract, and its deletion during development results in atresia of reproductive tract organs. Nearly all ovarian carcinomas express PAX8, and its knockdown results in apoptosis of ovarian cancer cells. To explore the role of PAX8 in these tissues, we purified the PAX8 protein complex from nonmalignant fallopian tube cells and high-grade serous ovarian carcinoma cell lines. We found that PAX8 was a member of a large chromatin remodeling complex and preferentially interacted with SOX17, another developmental transcription factor. Depleting either PAX8 or SOX17 from cancer cells altered the expression of factors involved in angiogenesis and functionally disrupted tubule and capillary formation in cell culture and mouse models. PAX8 and SOX17 in ovarian cancer cells promoted the secretion of angiogenic factors by suppressing the expression of SERPINE1, which encodes a proteinase inhibitor with antiangiogenic effects. The findings reveal a non-cell-autonomous function of these transcription factors in regulating angiogenesis in ovarian cancer.

Sex differences in microRNA expression in first and third trimester human placenta†.

Maternal and fetal pregnancy outcomes related to placental function vary based on fetal sex, which may be due to sexually dimorphic epigenetic regulation of RNA expression. We identified sexually dimorphic miRNA expression throughout gestation in human placentae. Next-generation sequencing identified miRNA expression profiles in first and third trimester uncomplicated pregnancies using tissue obtained at chorionic villous sampling (n = 113) and parturition (n = 47). Sequencing analysis identified 986 expressed mature miRNAs from female and male placentae at first and third trimester (baseMean>10). Of these, 11 sexually dimorphic (FDR < 0.05) miRNAs were identified in the first and 4 in the third trimester, all upregulated in females, including miR-361-5p, significant in both trimesters. Sex-specific analyses across gestation identified 677 differentially expressed (DE) miRNAs at FDR < 0.05 and baseMean>10, with 508 DE miRNAs in common between female-specific and male-specific analysis (269 upregulated in first trimester, 239 upregulated in third trimester). Of those, miR-4483 had the highest fold changes across gestation. There were 62.5% more female exclusive differences with fold change>2 across gestation than male exclusive (52 miRNAs vs 32 miRNAs), indicating miRNA expression across human gestation is sexually dimorphic. Pathway enrichment analysis identified significant pathways that were differentially regulated in first and third trimester as well as across gestation. This work provides the normative sex dimorphic miRNA atlas in first and third trimester, as well as the sex-independent and sex-specific placenta miRNA atlas across gestation, which may be used to identify biomarkers of placental function and direct functional studies investigating placental sex differences.